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Abstract. 

We use a novel parameterization of the flowing Hamiltonian to show that the flow 
equations based on continuous unitary transformations, as proposed by Wegner, can 
be implemented through a nonlinear partial differential equation involving one flow 
parameter and two system specific auxiliary variables. The implementation is non- 
perturbative as the partial differential equation involves a systematic expansion in 
fluctuations, controlled by the size of the system, rather than the coupling constant. 
The method is applied to the Lipkin model to construct a mapping which maps the 
non-interacting spectrum onto the interacting spectrum to a very high accuracy. This 
function is universal in the sense that the full spectrum for any (large) number of 
particles can be obtained from it. In a similar way expectation values for a large class 
of operators can be obtained, which also makes it possible to probe the stucture of the 
eigenstates. 
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1. Introduction and background 

The renormalization group and associated flow equations ^ is one of the cornerstones 
of modern physics and one of the very few potentially non-perturbative techniques 
available to us. Not only does the renormalization group allow the construction of 
low-energy effective theories for interacting quantum systems, but it also allows us 
to probe the phase structure of such systems. It is probably for this reason that 
considerable interest was generated in the flow equations recently proposed by Wegner j2] 
and independently by Glazek and Wilson [3j as well as its demonstrated close relation 
to the renormalization group 1^. Further applications of the flow equations, based 
on successive inflnitesimal unitary transformations, to diverse quantum mechanical 
problems followed, which included the treatment of the electron-phonon coupling [2], 
the spin-boson Hamiltonian |3], the Hubbard model [HI, the Sine-Gordon model [0] 
and the Foldy- Wouthuysen transformation [7j . The Lipkin model has been particularly 
prominent among applications [HI El EI], probably because it is a many-body model that 
can be solved numerically quite easily, but which exhibits a non-trivial quantum phase 
transition ^l] and could therefore shed light on the practical implementation of these 
flow equations in a simple, and yet non-trivial setting. 

One difficulty encountered in the practical implementation of the renormalization 
group and flow equations is that the Hamiltonian does not preserve its form under 
the flow and that additional operators, not present in the original Hamiltonian, are 
generated PlEj- The generic situation therefore yields an inflnite set of nonlinear coupled 
differential equations for the coupling constants, and some form of truncation is required. 
One way of doing this is through a perturbative implementation of the renormalization 
group equations, which corresponds to a systematic expansion in the coupling constant. 
Of course this procedure becomes invalid when flow to strong coupling occurs and non- 
perturbative aspects of the problem at hand cannot be probed in this way. In the 
context of the Wegner flow equation this issue has been dealt with in all the applications 
mentioned above by parameterizing the flowing Hamiltonian and then attempting to 
close on a flnite set of equations through some approximation. This procedure failed 
completely in situations where the system exhibited non-perturbative features, such as 
the quantum phase transition in the Lipkin model [S^, where the coupling simply flows 
to inflnity. A more detailed discussion on the role of truncation in the accuracy of the 
flow equation can be found in |12l. 

In a previous paper one of the present authors suggested a self-consistent 
implementation of the Wegner flow equation which yielded good results in both phases 
of the Lipkin model ^3] ■ The central idea of this approach was a systematic expansion of 
the flow equation in the fluctuations around the ground state expectation value, which is 
controlled by the inverse of the number of particles in the system, instead of the coupling 
constant. This resulted in an effective non-perturbative implementation of the flow 
equation, explaining its success in describing both phases of the Lipkin model. However, 
the approach as outlined in ^Sj made use of speciflc features of the Lipkin model and 
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was only successful in describing the low-energy behaviour, while it fails in producing 
the excitation energies of higher states correctly. Here we describe a generalization of 
this approach which can be applied to any finite dimensional Hamiltonian (see the next 
section for further elaboration) and which succeeds in producing the full spectrum and 
expectation values correctly. 

To set the scene and the notation we start with a brief review of the Wegner flow 
equation. The central idea behind the flow equation is to diagonalize the Hamiltonian 
with a family of successive infinitesimal unitary transformations, U{i), parameterized by 
a flow parameter i and generated by some appropriate anti-hermitian operator r]{i), i.e., 
'^^J'P = —U{i)ri(£). The resulting equation that determines the flow of the transformed 
Hamiltonian H{£) = W{£)H{0)U{i) is 

^ = W),H{i)]. (1) 

In Wegner 's original formulation the generator ?7(£) was chosen as the commutator 
of the diagonal part of the flowing Hamiltonian H{i) in some basis with the Hamiltonian 
itself, i.e., r]{i) = [diag(if(£)), H{i)] and it was shown that the final Hamiltonian H{oo) 
will be diagonal in the chosen basis . Here we allow a slightly more general formulation 
of the flow equation and consider a generator T]{i) = [Hq, H{i)] where it is assumed that 
the eigenvalues and eigenstates of Hq are known (note that Ho is i independent). It 
is easy to prove that the final Hamiltonian H{oo) is diagonal in the basis of Hq if the 
spectrum of Hq is non-degenerate, else it will assume a block diagonal structure. The 
proof rests on the following observation: 

±i,{H{i)-HQf = -2tTi[HQ,H{i)nHQ,H{i)]) < 0. (2) 

Thus tr(if(£) — Hq)^ is a monotonically decreasing function of £ which is bounded from 
below by zero. This implies that the derivative must vanish in the £ — > oo limit, and since 
the right-hand side of (j21) is simply the trace norm of [Hq, H{i)] this in turn implies that 
[Hq, H{oo)] = 0. Furthermore it can be shown that the eigenvalues of H = H{0), 
as they appear on the diagonal of H{oo), have the same ordering as the eigenvalues of 

Hq. 

In section 2 we describe the general parameterization on which our approach is 
based and outline the general derivation of the flow equation as an expansion in terms of 
the fluctuations. We briefly discuss the general properties of the resulting flow equation, 
which now takes the form of an ordinary (in contrast to operator valued) nonlinear 
partial differential equation, and describe how the spectrum and expectation values 
can be computed. In section 3 we apply this method to the Lipkin model where the 
calculations can be done simply and explicitly. In section 4 the results obtained by 
solving the resulting flow equation numerically are compared in detail with the exact 
results. The flow equation provides some interesting new perspectives on the phase 
transition and the general properties of the spectrum in the second phase. We conclude 
in section 5 with a summary, an outline of future applications and outstanding issues. 
The technicalities and proofs are given in three appendices. 
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2. A general parametrization of flowing operators 

We consider ttie flow of a real, hermitian operator H{i) which acts on a flnite dimensional 
Hilbert space Ti. Here we have in mind applications to systems with a flnite (but large) 
number of degrees of freedom such as interacting fermions or spins on a flnite sized 
lattice where one is in particular interested in studying the behaviour of the system as a 
function of size. For the moment we only consider systems in which there is no breaking 
of time reversal symmetry. Breaking of time reversal symmetry can be dealt with at 
the expense of a generalization of the expansion ()13|) and correspondingly the relations 
discussed below. 

The restriction to a flnite dimensional Hilbert space is made here since the expansion 
below is easily proven in this case (see Appendix A| ) . It is, however, by no means 



essential. If the Hilbert space is inflnite dimensional the expansion (jH)) still applies 
to bounded operators [T3] and one may study the flow of a bounded function of the 
Hamiltonian, instead of the Hamiltonian itself. Alternatively one can introduce a cut- 
ofi^ in some basis, e.g., a momentum cut-off and study the behaviour of the system as a 
function of the cut-off. 

As mentioned in the introduction the usual method of closing the flow equation 
consists of replacing H{i) by an approximate parameterized form, and then solving for 
the parameters as functions of i. This approach invariably fails at large coupling or when 
the system exhibits non-perturbative features as the approximate parameterization of 
H(£) breaks down. We abandon this approach in favour of a representation for H{i) 
which is both exact and very general. Let Hq and Hi be hermitian operators acting on 
7i which together form an irreducible set. What follows holds for any irreducible set, 
however, the case of two operators appear naturally in flow equations as the Hamiltonian 
can be written in the form H = Hq + Hi where the spectrum and eigenstates of Hq are 
known. Note that if Hq and Hi are reducible on H, but irreducible on a proper subspace 
of Ti, the problem can be restricted to this smaller subspace making the irreducibility 
of Ho and Hi a very natural requirement. In [Appendix A| it is shown that any operator 
acting on 7i can be written as a polynomial in Hq and Hi. In particular this holds for 
H{i) so we may write 

H{i) = r + rn, + r^H.Hj + r^^n^n^Hk + . . . , (3) 

where each F*-' " coeflicient is a function of £ and repeated indices indicate sums over 
and 1. Since H[tj is both hermitian and real it follows that the F's are invariant under 
the reversal of indices: 

rnin2...n, _ pn,nfe_i...ni > Q, G {0, 1}. (4) 

Now deflne Aifj = Hi — (Hi) as the fluctuation of Hi around the expectation value 
{Hi) = {(f) \Hi\ (f)) , where |0) is some arbitrary state. By setting Hi = AHi + {Hi) in 
equation (jS)) we obtain H{i) as an expansion in these fluctuations 

H{i) = f + VAHi + V^AHiAHj + V^^AHiAHjAHk + . . . , (5) 
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where 

f = r + r {H,) + {h,) {h,) + r^^ m {h,) {h^) + ... (6) 
f ^ = + (r*^' + P') (Hj) + [r^'' + P''' + P^') {Hj) (Hk) + ... (7) 

Y^j = + ^ ^ ^^^^ ^ _ _ _ 

and so forth. 

Note that T, when viewed as a function of i, (Hq) and (Hi), encodes information 
about all the expansion coefficients T^^'" appearing in equation Q- A natural strategy 
that presents itself is to set up an equation for F as a function of i, (Hq) and (Hi) in 
some domain T> of the {Hq)-{Hi) plane, determined by the properties of the operators 
Hq and Hi. In particular we require that this domain includes values of (-f^o) ranging 
from the largest to the smallest eigenvalues of Hq. For this purpose we require a family 
of states 10(A)) parameterized by a continuous set of variables X E V such that as A 
is varied over the domain V, (Hq) and (Hi) range continuously over the domain V. A 
very general set of states that meets these requirements are coherent states JH] • In this 
representation we may consider (Hq) and (Hi) as continuous variables with each f a 
function of {Hq), {Hi) and the following relations hold generally: 

P = -TTT-^ and F*^ = -—-4^——. (9) 

d{Hi) 2d{H,)d{H^) 

Coefficients with three or more indices cannot be written in this way, since, for example, 
r^°^ need not be equal to F*^^^. The general relationship between the true coefficients 
and derivatives of F are 

Qn+mY 



mm 



d {Hq)- d {H,)'' 

where the sum over p is over all the distinct ways of ordering n zeros and m ones. For 
example 

1^ = 1 (f + f 010 + r^oo) . (11) 

S\d{HQ)^d{Hi)' 3^ ^ ^ ^ 

Replacing f ooi by the left-hand side of the equation above is equivalent to approximating 
-pooi average of pooi^roio and F^oo general this amounts to approximating 

Yni...ni the average of the coefficients corresponding to all the distinct reorderings of 
ni, ...,ni. 

To set up an equation for T{{Hq) , {Hi) , i) we insert the expansion (0) into the flow 
equation (H)) with ri{i) = [Hq, H{i)] and take the expectation value with respect to the 
state 10(A)). The left-and right-hand sides become a systematic expansion in orders of 
the fluctuations (Aif„j . . . Aif„^): 
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^ + ^ (AHAH,) + ...= rP {[[AHo, AH,], AH,]) 

+ ff^'^' {[[AHo,AH,],AH,,AH,,]) 
+ PI^^F {[[AHo,AH,AH,^],AH,]) (12) 
+ f^'^T^^^'^ {[[AHo,AH,,AH,,],AH,,AH,,]) 
+ ... 

Note that writing Hq or AHq in the first position of the double commutator on the 
right-hand side of p2|) is a matter of taste. The expectation values appearing above are 
naturally functions of A and may be written as functions of (Hq) and (Hi) by inverting 
the equations (Hi) = (0(A) \Hi\ 0(A)) to obtain X{{Hq) , (Hi)). With an appropriately 
chosen state, such as a coherent state, the higher orders in the fluctuation can often 
be neglected, as the expansion is controlled by the inverse of the number of degrees 
of freedom (see Appendix C for an explicit example). Upon doing this we note that 



replacing the r^-^'^'s with more than two indices by a derivative will introduce corrections 
in fll2j) of an order higher than the terms already listed, or, on the level of the operator 
expansion, corrections higher than second order in the fluctuations. Working to second 
order in the fluctuations we can therefore safely replace F*-' " by the derivatives of a 
function / {{Hq) , {Hi) , i) = T{{Ho) , (H^) , £) and write for H{e): 

H{i) = f + f'AHo + f'AHi + i/ii {AHo, AH^} + ^fAH', + \rAHl + • • • 

= — —r (Distinct orderings of i AHq's and j AH[s) (13) 

ij=o ^^ + -^)- 

where {■, ■} denotes the anti-commutator and 

d {H,yd{HiY 

This turns the flow equation (fT^ into a nonlinear partial differential equation for 
/ {{Hq) , {Hi) correct up to the order shown in (fT^ . The choice of coherent state, 
the corresponding calculation of the fluctuations appearing in equation (fT^ and the 
identification of the parameter controlling the expansion are problem specific. Therefore 
we focus for the rest of this paper on a specific example, namely the Lipkin model, where 
all the calculations can be done explicitly. We would, however, like to emphasize that 
the results obtained here for the Lipkin model go well beyond any results previously 
obtained for the Lipkin or any other model within the context of the Wegner flow 
equations. 

Before embarking on the discussion of the Lipkin model, there are a number of 
general statements that can be made about the flow equation (|T2|l and the behaviour 
of f {{Hq) , {Hi) ,£). The first property to be noted is that since H{oo) is diagonal 
in the eigenbasis of Hq it should become a function of only Hq, provided that the 
spectrum of Hq is non-degenerate. This is reflected in the behaviour of / by the fact 
that / ((-f^o) ) {Hi) , i = oo) should be a function of {Hq) only. This is indeed borne out 
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to high accuracy in our later numerical investigations. This function in turn provides us 
with the functional dependence of H{oo) on Hq. Keeping in mind the unitary connection 
between H{oo) and H this enables us to compute the eigenvalues of H straightforwardly 
by inserting the supposedly known eigenvalues of Hq into the function / {{Hq) , oo). 

A second point to note is that the considerations above apply to the flow of an 
arbitrary real hermitian operator. It is easily verified that the transformed operator 
P{i) = W{e)P{0)U{i) satisfies the flow equation: 

^ = W),P{£)] (15) 

where the same choice of •)]{£) as for the Hamiltonian has to be made, i.e., •)]{£) = 
[Hq, H{i)]. The expansion (fT^ can be made for both operators H{i) and P{i). Denoting 
the corresponding function for P{£) by g {{Hq) , {Hi) ,£) the flow equation (fT3j) turns 
into a linear partial differential equation for g {{Hq) , {Hi) ,i) containing the function 
/ {{Hq) , {Hi) , £), which is determined by the nonlinear partial differential equation (jT^j). 
The expectation value of P = -P(O) in an eigenstate \En) of H can be expressed as 

(E„|P|E„) = {Er,\U{i)U\i)PU{i)U^{i)\E:,) = {E^,i\P{i)\Er„i) (16) 

where \En,i) = U''{i) \En). In the limit i ^ oo the states \En,oo) are simply the 
eigenstates of Hq, which are supposedly known. In this way the computation of the 
expectation value (P„ \P\ En) can be translated into the calculation of expectation 
values of the operator P(oo), obtained by solving the flow equation ()15|l . in the known 
eigenstates of Hq. 

3. The Lipkin Model 

The Lipkin model describes fermions distributed over two f2-fold degenerate levels 
separated by an energy of ^. The interaction V introduces scattering of pairs between 
shells. Labeling the two levels by a = ±1, the Hamiltonian reads 

H =^Y1 ^"-l^^P'^ ^ S (^laO'lf,a(^p',-aap,-a (17) 

where the sum over p,p' runs over the level degeneracy 1 . . .Q. A spin representation 
for H can be found by introducing the su(2) generators 

Jz = -'^cral^a(^p,a and J± = ^ (Ta],_±^ap,^„ (18) 

a,p p 

which satisfy [J^, J±] = ±J± and [J+, J_] = 27^. We divide H hj ^ and define the 
dimensionless coupling constant P = 2jV/^ to obtain 

H = J, + ^^{JI + J!) , (19) 

where all energies are now expressed in units of C,- Since [H, J^] = the Hamiltonian 
acts within irreducible representations of su(2), leading to a block diagonal structure 
of size 2j + 1. The low-lying states occur in the multiplet j = N/2. When (3 = the 
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ground state is simply \G) = a\ _^ . . . a^j^ _^ |0) which is written as \G) = \j = N/2, — j) 
in the spin basis. Non-zero values of (3 cause particle-hole excitations across the gap, 
and at /3 = ±1 the model shows a phase transition from an undeformed first phase to 
a deformed second phase. To distinguish the two phases we use the order parameter 
Q = l + {Jz) / i where (J^) is the expectation value of in the ground state. A discussion 
of this model and its features appear in ^T]. We use the shorthand \m) = \j,m) for 
the eigenstates of and denote the eigenstate of H with energy En by \En) where 
Eq ^ El < E2 .... 

3.1. Flow equations for the Lipkin Model 
Setting 

Ho = J. and Hi = ^ + J!) (20) 
the Hamiltonian becomes H = Hq + (3 Hi. We will consider flow equations of the form 

^^ = m,Hii)], HiO) = Ho + f3Hi (21) 
where the generator is chosen as in jSj, 

v{i) = [J.,H{i)] = [Ho,H{e)]. (22) 
The choice of in the generator has two important consequences: 

(i) Since Jz is non-degenerate and 77(00) = 0, H{oo) must be diagonal in the Jz basis. 

(ii) H{£) will remain band-diagonal, i.e. {n \H{i) \ m) = when \n — m\ ^ {0, 2}. 

The unitary transformation U {£) defined by 

^ = -UiiMi) and f/(0) = / (23) 

relates H{£) to H{0) in the usual way: 

H{i) = U\i)H{0)U{i). (24) 

The eigenvalues of H appear on the diagonal of if (00) in the same order as in Jz i.e. 
increasing from top to bottom. It follows that En = {—j + n \H{oo) \ — j + n) for each 
n = 0, . . . , 2j. The eigenstates of H{i) evolve as \En, i) = f^^(^) \En) during the flow, and 
since En = {En \H\ En) = {En, 00 |if(oo)| En, 00) we conclude that \En, 00) = \—j + n). 

3.2. Representing the flowing operators 

As suggested by the notation, we choose to represent flowing operators in terms of Hq 
and Hi as described in section [21 However {Hq, Hi} is not an irreducible set, since both 
leave the subspace spanned by {\n) \n even} invariant. As pointed out earlier this is 
easily remedied by restricting ourselves to this subspace. However, it turns out that in 
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Figure 1. The domain of f{{Ho) , (Hi) ,£). 



this case an expansion of the form Q can, as a consequence of the particular choice of 
the generator, even be found in the reducible case. Indeed, we may write H{i) as 

H{i) = Y,d'iW + Xl^^W^o^i^o- (25) 

i=0 j=0 

where the diS and rij's are the diagonal and off-diagonal coefficients respectively. This 
can be verified by inspection of the flow equation which shows that no new terms are 
generated. We will not use this model dependent parameterization in what follows, 
and only present it as proof that such a representation exists, demonstrating the more 
general theorem in this particular model. This result extends to f/(£) by equation ((221), 
and then to any flowing operator P{t) = W {£) P {0)U {£) , provided that P(0) can be 
expressed in terms of Hq and Hi . 

We will calculate the averages of Hq and Hi with respect to the coherent state 

\z)=N-^exp{zJ+)\-j), (26) 

where z = rexp{i9) G C and = 1 + r^. Using the method described in [Appendix B| 
it is found that 

,rr^ .fr^-l\ ^ ,rrx (2^ - l) COs(2^) 

W=j(^j and {Hi)= '' (/^^2)2 • (27) 

The possible values of {(Hq) , (Hi)) are shown in figure Q With these definitions in 
place we conclude that the flowing Lipkin Hamiltonian may be written as 

H{i) = f + fAHo + f'AHi + i/ii {AHo, AHi} + \rAHl + \rAHl + ■ ■ • (28) 

with Hq and Hi respectively the diagonal and off-diagonal parts of the original 
Hamiltonian. By definition AHi = Hi — {z\Hi \ z) and so / is defined on the domain 
pictured in figure ^ The initial condition becomes 

f{{H,) , {Hi) ,i = 0) = {Ho)+P (Hi) . (29) 
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Due to the coherent nature of the state (|26|) one might expect that terms corresponding 
to high order fluctuations in (PHj) will contribute less significantly to {z \H{i) \ z) than 
the scalar term /. This statement can be made precise as follows. Since the flowing 
Hamiltonian is extensive / should be proportional to j, as is the case with {Hq) and 
{Hi). To keep track of the orders of j we introduce the scaleless variables x = {Hq) /j, 
y = (Hi) / j and / = f / j. When taking the inner-product with respect to 1^;) on both 
sides of ()28|) the linear terms fall away and we obtain 



{Hm = f+ \r {{^Ho, AHi}) + \r {^hd + \r {^hd + .... m 

Each term is, up to a constant factor, of the form 

(Prod(^ + k)) = (Prod(z + k)) (31) 

where PTod{i + k) denotes some arbitrary product of i + k fluctuations. Using the results 
of Appendix C| we see that such a term is at most of order 0{j^^^^^''^j^^^^''^^'^^) in j. The 



leading order term corresponds to i + A; = 0, i.e. the scalar term /. We conclude that 
/ is the leading order contribution to {z \H{i) \ z), expressed not as a function of z and 
z*, but rather of the averages {Hq) and (-f^i). 

3. 3. The flow equation in the j ^ oo limit. 

We consider the flow of two operators, the Hamiltonian and another arbitrary operator 
P which flows as P{i) = U{i)''PU{i). Since H{i) determines U{i) one would expect a 
one-way couphng between the equations. It is assumed that P is a hermitian operator 
constructed in terms of Hq and Hi. Furthermore {z \P\ z) must be a rational function 
of j, which ensures that when taking derivatives of g{x,y,i) no additional factors of j 
are generated. First we summarize the equations concerned 

H{i) = f + ^AHq + f'AHi + i/ii {AHq, AHi} + ^fAH^, + ^^AHf + ■ ■ ■ ,(32) 
P(£) = g + g'^AHo + g^'AHi + ^g'' {AHq, AHi} + \9^'AHl + \g''AHl + ■ ■ ■ , (33) 
'^^ = [[HQ,Hm,Hm, H{0) = Hq + PHi, (34) 

^ = [[HQ,H{i)],P{i)], P(0) = P, (35) 

where /{{Hq) , {Hi) , 0) = {Hq) + P {Hi) and gi{HQ) , {Hi) , 0) is just {z \P\ z) up to 
leading order in j. Next we substitute the expansion of H{i) into the flow equation ()34j) 
and take the expectation value on both sides with respect to the state \z). Arguing as 
before we identify the leading order term on the left as being df / di. A general term on 
the right is of the form 

By transforming to scaleless variables and using the result of Appendix C it is seen to 



be at most of order 0{j^ *+L(i+*)/2j^^ where t = i + n + k + l. Since the t = and t = 1 
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terms are zero (they involve commutators of scalars) the leading order contributions 
come from the t = 2,3 terms. These are exactly the terms found by considering the 
expansion of H{£) up to second order in the fluctuations. The expectation values of the 
double commutators may be calculated and expressed as functions of x and y using the 
method outlined in [Appendix B It is found that, due to cancellations, only four of the 



potential fifteen terms make leading order contributions. The corresponding expectation 
values are, to leading order 

{[[AHo,AH,],AHo]) = g2iix,y)j (37) 

{[[AHo,AH,],AH^]) = 922{x,y)j (38) 

{[[AHo,{^Ho,AH,}],AH,]) = g,,{x,y)f (39) 

{[[AHo, AH^] , AHo] > = - (732(0;, y)f (40) 

where g2i = —4:y, (722 = 2x(l — x"^), g^2 = ~2 + 4x^ — 2x^ + 8?/^. When keeping only the 
leading order terms on both sides the flow equation becomes 

% = 92.~r~r + 922 ~r~r + ^^?32/^v~°^ - ^^732/°v~^°, (41) 

where superscripts denote derivatives to the rescaled averages x and y. For further 
discussion it is convenient to transform to the variables (p = r^,g = cos(26')) G 
[0, 00) X [—1, 1], which are related to {x,y) by 

a: = P^ and y = ^^ + 0(i). (42) 

Note that the domain of these variables is a strip, which significantly simplifies the 
numerical solution. We note that the arguments above can be applied, completely 
unchanged, to the flow equation of P{t) as well. We obtain, now for both H{t) and 
P(£), the coupled set 

— = -2(1 + P) + (1 - 9 (43) 



^ Y dq dp ^ dq dqdp j 

di '^i^ ~^ P) \^dqdp~^^^ ^ dq dqdp I ' 



(44) 



Note that in contrast to the equation for /, the equation for g is a linear equation that 
can be solved once / has been obtained from ()43|) . 

In section l!m it was mentioned that H{i) retains its band diagonal structure during 
flow, which meant that Hi appeared only linearly in the representation of H{i). This 
implies that f{{Ho) , (Hi) , £) should be linear in (Hi), or, in the new variables, linear in 
q. When the form f{p,q,i) = no{p,i) + qni{p,i) is substituted into (jlH)) this is indeed 
seen to be the case, and we obtain a remarkably simple set of coupled PDE's for uq and 
ni: 

-di=-^^^+P^ "^9^ 

drii , ,9 dnn , , 

_l = _2(i+,)V^. (45) 
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The initial conditions are 

„„(p,0) = .= f^ and „,(p.0) = * = ^. (46) 

p+l Q [P + V 

The matrix elements of the Lipkin Hamiltonian possess the symmetry {m \H\m) = 
— {—m\H\ —m), which implies that the spectrum is anti-symmetric around Ej = 0. 
This symmetry is respected by the flow equation and manifests itself through an 
invariance of equations under the substitutions Uq —no, rii rii and p p~^. 
For the solutions this implies that nQ{p,i) = —nQ{p~^,i) and ni{p,i) = ni{p~^,(^), and 
we may thus restrict ourselves to the interval [0, 1] in the p dimension. This will be seen 
to correspond to the negative half of the spectrum, from which the entire spectrum can 
easily be obtained. Apart from the initial conditions we impose the boundary conditions 
no{l,i) = 0, as required by symmetry considerations, and ni(0,£) = 0, since at r = 
the 6 dependency should vanish. The latter is also consistent with and, in fact, follows 
from the flow equation and initial value at p = 0, since dni/di oc ni. A naive application 
of the same argument to no seems to suggest that no flow occurs for no at p = and 
that rzo(0, = —1 for all i. This conclusion is, however, incorrect since a more careful 
analysis shows that in the second phase drii/dp develops a square root singularity at 
p = 0, which allows no(0,£) to flow away from —1. Numerically this can be handled 
easily by solving for m(p, £) = n\{j), £), instead of rii. 



3.4- Finding the spectrum 

From / = j/ it is possible to obtain the entire spectrum of H. For this discussion it 
is convenient to use the original {{Hq) , (Hi)) variables. Recall that in the / — > oo limit 
H{£) flows toward a diagonal form, furthermore the eigenvalues appear on the diagonal 
in the same order as in Hq, i.e. increasing from top to bottom. This implies that the 
^th eigenvalue of H is given by 

En = {-J + n \H{£ = oo)\ -j + n), (47) 

where Eq corresponds to the ground state energy. As H{£) flows towards a diagonal 
form the terms of expansion (jH^ containing Hi will disappear and eventually H{oo) 
and / will become functions of Hq and (Hq) only. The eigenvalues of H are given by 

En = fi{Ho) = -J + n, 0, £ = oo), (48) 

where (Hi) has been set to zero. This can be understood in two ways. Looking at 
equations (jHl) and © we see that the functional dependence of H{£^) on Hq is the same 
as that of T (now called /) on (Hq). Since taking the expectation value of H{oo) 
with respect to |— j + n) is equivalent to substituting —j + n for Hq, the result follows. 
Alternatively, consider equation (fT!?j) and note that setting (Hq) = —j + n makes the 
(n + 1)^* diagonal element of AHq zero. Since Hi does not appear we see that when 
taking the inner product with \—j + n) only the scaler term, /, will survive. The results 
of this procedure appear in section 14.11 
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3.5. Calculating expectation values 

Next we return to the arbitrary operator P introduced earlier, and show how its 
expectation values with respect to the eigenstates of H may be calculated. As examples 
we consider the first and second moments of Hq = in the ground state. The 
former is simply related to the order parameter fl discussed in section El Suppose 
P is a hermitian operator having the properties mentioned in section 13.31 The aim 
is to calculate {En \P\ En) where \En) is the (unknown) eigenstate of H = Hq + PHi 
corresponding to the energy En- We can formulate this calculation in terms of a flow 
equation by noting that 

{En \P\En) = {En \U {e)U\e) PU {e)U\i)\ En) = {En,i\P{i)\EnJ) (49) 

where \EnJ) = W{i)\En), P{e) = W{i)PU{i) and U{i) is the unitary operator 
associated with the flow of H{£). This equation holds for all i, and particularly in 
the i ^ oo limit. Since \En, oo) = + n) it follows that 

{En \P\ En) = {-J + n \P{oo)\ -j+n), (50) 

and we conclude that the expectation value of P in the n^^ eigenstate of H is simply 
the n^^ diagonal element of P(oo) in the Hq basis. Furthermore {z\P{i)\z) has the 
character of a generating function in the sense that 
Qn+m ^z\P{e)\z) 



(x{-j+n\P{i)\-j+m), (51) 

•=o 

thus knowing g, which is just {z \P{i) \ z) up to leading order in j, is sufficient to obtain all 
the matrix elements of P{i) with high accuracy. However, unless an analytic solution 
for g is known, we are limited to numerical calculations for the low lying states. In 
particular, the ground state expectation value is found by setting z = 0. 

Consider the case where P = Hq/J = Jz/j and g{x,y,i = 0) = x. Here P{£) will 
flow from a diagonal form to one containing high powers of Hi. Having no knowledge 
about the form of g we must solve the flow equation directly as a linear PDE in 3 
variables. For numerical implementation we again change to the variables p, q and 
solve on the strip p G [0, 1], q G [—1, 1]. As before we impose the boundary condition 
(||)p=Q = 0. It turns out that the solution for p << 1 is virtually independent of the 
boundary condition imposed at p = 1 (this is in fact also true for / or uq and ni) and 
in this case we simply allow the solution to evolve freely at p = 1. The results of this 
calculation is shown in section 14.21 



3.6. Probing the structure of eigenstates 

The ability to calculate expectation values for the class of operators discussed above 
allows us to probe the eigenstates in a variety of ways. We will consider diagonal 
operators of the form 0{h) = Yln=-j ^(^/j) 1^) (^1 where h is a. smooth, j-independent 
function defined on [—1, 1]. Since Hq is non-degenerate its powers span the space of 
diagonal operators, and so 0{h) can always be expressed as a function of Hq. In order 
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to apply the flow equation to 0{h) we must find the initial condition for equation (j44p . 
which is given by 

g{x,y,i = 0) = hm {z\0{h)\z) . (52) 

The expansion coefficients for the coherent state in the Jz basis are given by ^\ 

2j ^ 



{n\z) = (1 + r')-^ .] z-+^ n = -j,...,j (53) 

\n + jj 

which, together with a change of variables from r to x = (r^ — l)/(r^ + 1), lead to 

This is recognized as the Bernstein polynomial jlHj for h{x) of degree 2j on the interval 
[—1,1]. It is well known that in the j ^ oo limit this polynomial approximation con- 
verges uniformly to h{x), from which we conclude that the boundary condition of the 
flow equation is simply g{x,y,i = 0) = h{x). 



Suppose \Em) = J2n=-j |^) is the eigenstate under consideration, in which case 
{E^\0{h)\E^)= Y,h(-\\an?= 'f^h(-\\{n\E^)\\ (55) 

n=-j V-'/ n=-j ^-^ / 

This leads to the interpretation of the expectation value as an average of the expansion 
coefficients | {n\Em) p, weighted by the function h{x). It may appear that choosing 
0{h) to be the projection operator onto some basis state \n) provides the most direct 
means of calculating the contributions of individual states. However the projection 
operator does not fall within the class of operators for which equation (jl^ was derived. 
This follows from the fact that when taking derivatives of {z\n) {n\z) new factors of j 
are generated, which is equivalent to the observation that no continuous function h{x) 
exists such that 0{h{x)) = \n) {n\ for all j. Instead we choose h{x, x) = exp(— 7(x — x)^) 
where x G [—1,1] and 7 >> 1. This weight function focuses on the contribution of those 
basis states \n) for which n/j lies in a narrow region centered around x. Considered as 
functions of x, it is expected that {0{h{x,x))) would approximate limj^oo \c(jxf up to 
a constant factor, provided that the latter varies slowly on the scale of I/7 in the region 
of X. Clearly 7 controls the accuracy of this method, and also determines the scale on 
which the structure of the eigenstates can be resolved. 

4. Numerical Results 

4.1. The flow ofH{l) 

In section it was shown the (n + 1)*^* eigenvalue of H is given by = / (— j + n), 
where we suppress the last two arguments and take {Hi) = and £ = cxd in what follows. 
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Figure 2. Results obtained from the numerical solution of the flow equations, (a) 
The full spectrum. For clarity the exact result is only displayed for every fifteenth 
eigenvalue, (b) Relative error in the first five states as a function of j for /? = 0.5. (c) 
The order parameter = 1 + (Hq) /j (in the ground state) as a function of /3. (d) 
(i?o)^ /i^ in the ground state as a function of /?. 



En^jf{x^-1 + n/j) = jno p 



= jno (x) , 



(56) 



Note that / may be interpreted as a mapping of the non-interacting spectrum onto the 
interacting spectrum of H at finite /3. In terms of the rescaled variables this becomes 

1-x 

where we use ni(p,£ — oo) = 0. As j ^ oo the difference between values of x 
corresponding to successive states become zero, and we may consider x e [—1, 1] as 
a continuous label for the eigenstates through the assosiation x <-> \Ej(^i+x)}- 

Another property of interest is the gap between successive levels, which, for large 
j, is given by 



En+l — En — j 



no X = -1 + 



n+1 



- no X 



-1 + 



n 



dno 
dx 



(57) 



-l+n/j 
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Figure 3. Energy levels and gaps as functions of the coupling constant, (a) The ground 
state energy, (b) The gap between the second excited state and ground state, (c) The 
gap between the fortieth excited state and ground state, (d) Gap between the first 
excited state and ground state. 



Figure |H (a) shows the numerical solution for no{x) at £ = 6 and for three different 
values of (3. The inset is the corresponding graph for log(ni). For these values of the 
flow parameter rii, which represents the off-diagonal part of H{i), is already of order 
10^^'^ and may be disregarded completely. As a first test we calculate the spectrum for 
j = 1000 from direct diagonalization and plot the pairs (x = — 1 + n / j , E^'^'^^ / j) for 
n = 0, . . . , j as dots. There is an excellent correspondence for all states and values of 
the coupling constant as reflected in figure |21 (b) which shows the relative error in the 
first five energy eigenvalues as a function of j for (3 = 0.5. The lines fall almost exactly 
on one another, so no legend is given. The log-log inset shows a set of straight lines 
with gradients equal to one, clearly illustrating the 1/j behaviour of the errors. The 
same results hold in the second phase, although here the eigenvalues corresponding to 
the even states n = 0, 2, 4, . . . tend to be much more accurate than those for odd states. 
We return to this point shortly. 

Next we investigate how the eigenvalues depend on the coupling [3. Figure El (a) 
shows the ground state energy as a function of (3 for j = 1000 together with the exact 



Non-perturbative flow equations 



17 



result. The relative error is about 0.08%. As one would expect from the global spectrum 
in figure |21 the gaps are produced to a good accuracy. As examples we display the gaps 
E2 — Eq and £^40 — Eq in figures El (b) and (c). The absolute error is shown as an inset. 

One quantity that is not correctly reproduced is the gap between the ground state 
and first excited state. Figure El (d) displays this quantity, which turns out to be correct 
only in the first phase. Exact diagonalization shows that this gap vanishes in the second 
phase as 1/j, while the fiow equation produces a gap growing linearly with (3. The 
origin of this problem is, however, well understood and relates to the observation made 
above regarding the accuracy of the even states n = 0, 2, 4, . . . versus the odd states in 
the second phase. The reason why the fiow equations, at least as implemented here, fail 
to reproduce this quantity correctly relates to the lack of irreducibilty of {Ho,Hi} on 
the full space (see section IT2|) . In principle one should implement the fiow equations 
on each invariant subspace {\n) \n even} and {|n) \n odd} seperately, in contrast to 
what has been done here. This will produce two functions, mapping the non-interacting 
spectrum in each subspace onto the corresponding interacting spectrum. In this way the 
correct features of both sets of excitations can be reproduced. This demonstrates that 
care should be taken to ensure that the underlying symmetries and associated selection 
rules of the Hamiltonian are properly taken into account in the fiow equation. 

The fiow equation provides an interesting new perspective on the quantum phase 
transition which we now discuss. We mention that the graphs in figured which illustrate 
properties of the Lipkin model not directly related to fiow equations, were obtained from 
the direct diagonalization of H for j = 1000. Where applicable the fiow equation gives 
the same results. As is well known [19^ a quantum phase transition is related to a large 
number of avoided level crossings. This is clearly displayed in figuresEl(b) and (c) where 
the approximate vanishing of the gap at the critical value (3 = 1 can clearly be seen. 
From equation (jK7|l we note that a vanishing gap corresponds to a vanishing derivative 
of no(x). From this point of view, avoided level crossings manifest themselves as a 
vanishing derivative of no{x), the function that maps the non- interacting spectrum onto 
the interacting spectrum. From figure |21 (a) it is clear that such a vanishing derivative 
never occurs when (3 < 1. For values of /3 > 1, however, one always finds a point where 
the curve no{x) forms a plateau and the derivative vanishes approximately as shown 
in figure El (a). In fact at /3 = 1 this happens right at x = —1, i.e., at the ground 
state energy so that at this value of /3 the gap between the lowest excited states and 
the ground state vanishes, consistent with our earlier observation. At larger values of 
/3 this point shifts up in the spectrum as is clear from figure |2l (a). Note from figure 
121 (a) that at this point no = —1 for all values of f3. We remark that the approximate 
vanishing of the derivative is probably due to numerics. If the derivative does vanish 
exactly at some point, which we beheve to be the case, it will be very hard to detect it 
numerically as the fiow equation is of course solved on a grid with finite spacing. Note 
that a vanishing derivative of no(x) does not imply a vanishing gap at finite j (number of 
particles) as no{x) is then only evaluated at a finite number of discrete points, excluding 
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to all probability the point where the derivative of no(x) vanishes. It must be kept in 
mind that no(x) really reflects the situation in the thermodynamic limit and therefore 
it displays a point where the derivative vanishes. 

Returning to the flow equation ()45p we make the interesting observation that at 
the point where the derivative of uq^x) vanishes, ni(x) is not forced to flow to zero, 
but may in fact attain a non-trivial fixed point value. This is demonstrated in the inset 
of figure |21 (a) which shows log(ni(a;)) as a function of x. One clearly sees a sharp 
peak at the point where the derivative of no{x) vanishes approximately. This peak 
only occurs for f3 > 1 and moves to the right as f3 is increased. This is also consistent 
with the known result that an off-diagonal element rriij of H{i) decays roughly as 
exp(— — Ej_i)^i). This suggests a connection between quantum phase transitions, 
the corresponding disappearance of an energy scale (gap) [19] in the thermodynamic 
limit and the absence of decoupling in the Hamiltonian, also in the thermodynamic 
limit. 

Upon further analysis of the flow equation one observes that the point where the 
derivative of uqIx) vanishes, must necessarily be a point of inflection. The reason for 
this is that the derivative of no{x) cannot change its sign as the flow equation for ni{x) 
will become unstable upon which ni{x) will grow exponentially, something that cannot 
happen as explained in section^ (see equation |21) . This explains why we only observe a 
plateau in figure El (a) and no more drastic behaviour. Note, however, from figure ID (a) 
that, as can happen at a point of inflection, the second order derivative of no{x) does 
not exist, but changes sign. The fact that the dervative of no{x) cannot change sign 
also explains the remark often made that the ordering of the eigenvalues is the same in 
the interacting and non-interacting systems. 

Against the background of the discussion above, the presence of a point of inflection 
at higher energies for values of j3 > 1 presents the interesting interpretation of a 
'quantum phase transition' at higher energies. To illustrate that this is indeed a 
sensible interpretation we show the expectation value 1 + (J^) /j at /? = 0.5,1.5,2.5 
as a function of evaluated in the excited states, in figure lU (b). One clearly 

sees a change in the behaviour at the point of inflection. Indeed, after this point the 
behaviour of this quantity is linear, as is to be expected in the first phase, so that one 
can think of this point as a transition from the second phase, present at lower energies, 
to the first phase at still higher energies. This interpretation is summarized in a 'phase 
diagram' in figure ID (c) which shows the (3 dependancy of a subset of the negative, 
even eigenvalues. For (3 < 1 the eigenvalues are confined between and —j. As (5 
increases first the groundstate (shown in bold) and then the excited states begin to 
cross the E = —j phase boundary until eventually, in the large coupling limit, only 
the \Ej = 0) state retains its first phase character. For finite j successive eigenvalues 
show avoided level crossings on the phase boundary. In the thermodynamic limit one 
would, however, expect that successive eigenvalues will coalesce as they cross the phase 
boundary; signalling a vanishing gap. Keeping in mind the symmetry En = —E2j-n, 
these remarks can easily be adapted for the positive half of the spectrum as well. 
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Figure 4. (a) Derivative of no{x) as a function of x for /3 ~ 1,2.5. (b) 
1 + (£^j(i+a;) \ Jz\ Ejf^i^x)) 1 3 as a function of a; for /3 = 0.5,1.5,2.5 and j = 1000. 
(c) The negative half of the spectrum as a function of /3 for j ~ 1000. Only every 
twenty-fifth even eigenvalue is shown. 



The order parameter 

Next we present results for the expectation values, computed according to the method 
described in section 13.51 Figure |21 (c) shows the order parameter i7 = 1 + (i^o) 1 3 
as a function of /3 together with the exact values for j = 1000. Again we find 
excellent agreement, even at the phase transition, with absolute errors in the range 
of 5 X 10~^. The second moment of is found by choosing P = Jl/j"^-, for which 
giiHo) , (Hi) ,i = 0) = (Hq)'^ / j"^. Proceeding as before we obtain the expectation value 
of P as a function of (3, which appears, together with the exact result for j = 1000, in 
figure El (d). Note that the computation of the expectation values of different operators 
entails solving the same partial differential equation, but with different initial conditions. 
In this sense the flow equation only knows about the operator being computed through 
the initial conditions. 
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X 

Figure 5. The expectation values of 0{h{x, x')) at f3 — 1.3 for different x' . All results 
have been normalized to have the same maximum value. 

4.3. The structure of the eigenstates 

We follow the strategy discussed in section and consider the flow of the operator 
0{h{x,x)) for h{x,x) = exp(— 7(0; — x)^) and 7 = 500,2000. The resulting ground 
state expectation values for j3 = 1.3 and a range of x values appear in figure The 
exact values of \ {xj\Eo)\, calculated from direct diagonalization, are shown as functions 
of X. Since the ground state lies within the subspace spanned by even basis states 
we completely neglect the coefficients corresponding to odd states. As before the 
flow equations cannot distinguish these two subspaces on this level, and only provide 
information on the average contribution of a range of both even and odd basis states. 
As expected, larger values of 7 provide a more accurate probing of the eigenstate, and 
for 7 = 2000 we obtain a good approximation for the form of the absolute wave function 
at large j in the Jz basis. 

In the first phase a very small subset of basis states has nonzero overlap with the 
ground state. This implies that, as a function of x, limj^oo K^^jI-^'o)! will vary on a 
length scale which cannot easily be probed using this method. However, by considering 
0{h{x) = x) = Jz/ i and keeping in mind that {Eq \ Jz\ Eq) = —j \Eq) = we 
conclude from figure|21(c) and (d) that \Eq) is approximately equal to |— j) when (3 < 1. 

5. Summary and outlook 

We have shown that the fiow equations based on continuous unitary transformations 
can be cast in the form of a general nonlinear partial differential equation involving 
one flow and two system speciflc auxiliary variables. This approach is based on the 
observation that operators acting on finite dimensional Hilbert spaces, or bounded 
operators acting on infinite dimensional Hilbert spaces, may be written as polynomials 
in a set of irreducible operators, specifically two in the case of the fiow equation. 
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Making an expansion to second order in the fluctuations yields the aforementioned 
flow equation. It is shown how the spectrum and expectation values can be extracted 
from the solution of this equation. The method is applied to the Lipkin model where 
all the calculations can be done easily and explicitly. Excellent results for the full 
spectrum and expectation values are obtained. It was illustrated how properties of the 
eigenstates can be investigated by considering the flow of certain diagonal operators. 
Further information about the eigenstates, for example the relative phases of expansion 
coefficients, could possibly be obtained by extending this strategy to off-diagonal 
operators, or by considering the flow of U{i) itself. The flow equation also provides 
some interesting new insights into the nature of the quantum phase transition. 

A further possible extension of this approach is to study time dependent problems. 
Of course, once the eigenvalues and eigenfunctions have been obtained to a reliable 
accuracy one can, in principle, determine the time dependency of the system. We 
would, however, like to flnd a more economical way of computing the time dependency. 

Further possible applications to lattice models are currently being investigated. 
In particular one-dimensional lattice models should be amenable to this approach as 
they are in structure very similar to the Lipkin model. Indeed, the Lipkin model can 
be viewed as a lattice model with position dependent hopping amplitude and a linear 
potential. 
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Appendix A. Representations by irreducible sets 

In section|21we made use of the fact that any flowing operator may be expressed in terms 
of operators coming from an irreducible set which contains the identity. Now we show 
why this is the case, using a result by von Neumann. First we establish some notation. 
Let 7i be a flnite dimensional Hilbert space and denote by B{T-l) the set of all matrix 
operators acting on it. The commutant Ai' of a set C B{T-l) is deflned as the set of 
operators which commute with all the elements of Ai. The double commutant of Ai is 
deflned by At" = {At'}' . The following theorem [20] provides the key: 

The Double Commutant Theorem: If A is a subalgebra ofB{T-C) which contains 
the identity and is closed under hermitian conjugation then A = A" . 

Now suppose S C i3(7i) is an irreducible set of hermitian operators, one of which 
is the identity, and A the subalgebra of B(H) spanned by all the products of operators 
from S. By Schur's lemma S' = A' = {XI | A G C} and thus A" = B(H). It follows from 
the theorem above that A = BiH), i.e. any operator acting on B(TC) may be written in 
terms of operators coming from S. 
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Appendix B. Calculating expectation values with respect to coherent states 

Define l^;) = (1 + zz*)~^\z) = N~^z) where \z) = exp(zJ+) is the unnormalized 

coherent state. The differential operator representations for the su(2) generators are 
given by [T7] : 

W= (-J + \z) = I {z, d,) \z) (B.l) 

J+\z) = -^Jz) = J+{z,d,)\z) (B.2) 

J-\z)= (^2zj-z^-^^ \z) = J.{z,d,)\z). (B.3) 

We will always consider z and z* to be distinct variables, and that (^1 is a function of 
z* only. When calculating expectation values of the double commutators in section El 
one encounters expressions of the form 



{z\H„,H^,...Hr.Jz) = ^"'''"^ ,"V; (B-4) 



(^1 




■ ■ ■ Hnc, 




(^1 


\Z) 



where nj = or 1. These can easily be computed by replacing each operator by its 
differential representation to obtain 

{z\Hn^Hn2--Hn^z) = -^--H„^ {z,d^) . . . {z,d,) {z\z), (B.5) 

{z\z) 

where the orders of the operators have been reversed. As a result these expectation 
values can be expressed as some rational function of z = rexp{i6) and z* = rexp{—i6). 
Applying this method to (Hq) and (Hi) we obtain 



2 



{z\Ho\z)=j\^-^j and {z\H,\z) = (^^r^ip ■ (B-6) 

Using these equations to express r and 6 in terms of (Hi) we may consider expectation 
values of the form of equation ()B.4|) as functions of these averages. For example 

{z\[[AHo,AH^],{AHoy]\z) = 8Jxy, (B.7) 

where (Hq) = jx and {Hi) = jy. 



Appendix C. Scaling behaviour of fluctuations 

We wish to show that 

{z\AHn,AHn2...AHn,\z) j^^^ for V t G N whcrc G {0, 1} . (C.l) 

From direct calculation this is found to hold for t = 2, and we employ induction to 
obtain the general result. Assuming that this holds for all products oi k < t fluctua- 
tions, the induction step consists of adding either AHq or AHi to a general product 
M = AHn^AHn2---AHrn and then proving the result for {z \ MAHi \ z). 
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First we add an extra AHq. Using the results from the previous section we may 
write 

{z\MAHo\z)= - j{z\M\z) + {z\Mzd,\z) - {Ho){z\M\z) (C.2) 
= -2jzz*{l + zz*)-\z\M\z) ~ z{z\{d,M)\z) + zd,{z\M\z) (C.3) 

where ()B.6jl and {z\Mzdz\z) = z{dz{z\M\z) — {z\{dzM)\z)) were used. It should be 
remembered that M contains a z dependency through the average appearing in each 
AHi. Dividing by {z\z) leads to 

{z\MAHo\z) = -2jzz*{l + zz*)-^{z\M\z) - z{z\{dzM)\z) 

+ ^l + l,*y, dMM\z). (C.4) 

The last term can be rewritten using 

^^r^^dz{z\M\z) = dz{z\M\z) + (YT^^^I^I^)' (C.5) 
which is just the product rule, to obtain 

{z\MAHo\z) = zdz{z\M\z) - z{z\{dzM)\z). (C.6) 

Note the important cancellation of terms proportional to j{z\M\z). From the induction 
hypothesis, and the fact that {z\M\z) is polynomial in j, the first term in (j(I6j) will be 
of order j LtJ . The second term becomes 

t 

zJ2{^\AK,,...AH^^_^AH^^^^...AH^^z){dz{H^^)). (C.7) 

i=l 

Taking into account the expressions for (Hq) and (-f^i), we conclude that each term in 
equation ()C.7|) will have order j^^^^^, which reduces to j'-^-' in both the cases where 
t is odd and even. Thus, for any product M of t fluctuations we arrive at 

{z\MAHo\z) ^ j^'-^^, (C.8) 

which concludes the induction step. 

Exactly the same procedure is followed when adding a AHi, although more algebra 
is required as Hi{z,dz) now contains second order derivatives to z. The final result 
remains unchanged: 

{z\mAHi\z) ~ jL^J. (C.9) 

In this case these results make exact the general notion that relative fluctuations 
scale like powers of one over the system size. Finally we mention that when calculating 
the expectation values of a double commutator there is often a cancellation of leading 
order terms. This can be seen in equation ()B.7|) . where the sum of terms of order j'^ 
turns out to be of order j. 
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